Theory of the anomalous Hall effect from the Kubo formula and the Dirac equation 
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A model to treat the anomalous Hall effect is developed. Based on the Kubo formalism and on 
the Dirac equation, this model allows the simultaneous calculation of the skew-scattering and side- 
jump contributions to the anomalous Hall conductivity. The continuity and the consistency with the 
weak-relativistic limit described by the Pauli Hamiltonian is shown. For both approaches, Dirac and 
Pauli, the Feynman diagrams, which lead to the skew-scattering and the side-jump contributions, 
are underlined. In order to illustrate this method, we apply it to a particular case: a ferromagnetic 
bulk compound in the limit of weak-scattering and free-electrons approximation. Explicit expres- 
sions for the anomalous Hall conductivity for both skew-scattering and side-jump mechanisms are 
obtained. Within this model, the recently predicted "spin Hall effect" appears naturally. 



I. INTRODUCTION 

The Hall resistivity of magnetic materials, in addition 
to the normal part proportional to the magnetic field, 
contains a supplementary part proportional to the mag- 
netization, called the anomalous Hall resistivity 



p H = RoH + R S M, 



(1) 



where i?o and R$ are the normal and anomalous Hall 
coefficients respectively, H the magnetic field and M the 
magnetization. While the normal Hall effect results from 
the Lorenz force, the anomalous Hall effect is due to the 
spin-orbit coupling in the presence of spin polarization. 
Experimentally, the normal and anomalous parts can be 
extracted by measuring the Hall resistivity as a func- 
tion of the magnetic field. At high magnetic field, when 
the magnetic saturation is reached, we get a linear vari- 
ation of the Hall resistivity with a slope related to i?o 
and an extrapolated value at zero magnetic field related 
to Rs- The normal and anomalous Hall coefficients have 
been determined for a large number of bulk alloys. These 
studiestrQ reveal that the sign of Rs can change accord- 
ing to the alloy composition and that |i?sM| is generally 
larger than \RqH\ for typical values of the magnetic field. 

For different reasons, renewed attention to the anoma- 
lous Hall effect is observed quite recently. It is not only 
due to the increasing interest in spin-dependent trans- 
port phenomena but also because of some particular and 
interesting behaviors of the anomalous Hall resistivity 
obtained experimentally in granular alloysB, in magnetic 
filmsD and multilayersO. In addition, the anomalous Hall 
effect is increasingly used as a measurement tool to de- 
tect for .example magnetizationLl, dynamics of magnetic 
domainal!] or perpendicular anisotropylH Besides, a new 
effect closely related to the anomalous Hall effect, the 
"spin Hall effect" , has been recently predictec 

In the sixties, a number of theoretical worksl 
tempted to elucidate the physical mechanisms respon- 
sible for the anomalous Hall effect and to calculate an 
explicit expression forr-thp-|anomalous Hall resistivity. A 
series of controversiesc3~t3 arose from those pioneering 




works which were-soived through detailed calculationsEll 
and comparisonsEll. It is now acceptedna that two mech- 
anisms are responsible for the anomalous Hall effect: the 
skew-scattering proposed by Smito and the side-jump 
proposed by Bergenia. 
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FIG. 1. Schematic picture of the skew-scattering (a) and 
side-jump (b) mechanisms from a quantum point of view (0 
corresponds to spin up and to spin down) . The bold curves 
represent the anisotropic enhancement of the amplitude of the 
wave-packet due to spin-orbit coupling. 

An illustrative picture of these mechanisms is given in 
Fig. 0. Consider an incident plane-wave characterized by 
a wave- vector k which is scattered by a central potential 
due, for example, to impurity. In the presence of spin- 
orbit coupling, the amplitude of the wave-packet becomes 
anisotropic in the sense that it depends of the relative di- 
rections of the scattered and incident waves and of the 
spin. After a succession of scattering events, the average 
trajectory of the electron is deflected by a spin-dependent 
angle, which is typically of order 10~ 2 rad. This first 
mechanism, depicted by the diagram (a) in Fig. |], cor- 
responds to the skew-scattering. The second mechanism 
corresponds to a lateral displacement, S « I0~ n m, of 
the center of the wave-packet during the scattering, which 
is also spin-dependent. This mechanism, depicted by the 
diagram (b) in Fig. ^ corresponds to the side-jump. In 
both cases, due to the spin-orbit coupling, the effect is 
asymmetrical in respect to the spin state. The spin up 
and spin down currents are then different. In magnetic 
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materials, this leads to a non-zero spin current and to a 
transverse component in the charge current, which cor- 
responds to the anomalous Hall effect. 

The skew-scattering and the side-jump mechanisms 
give different contributions to the anomalous Hall resis- 
tivity. For bulk material, it has been shown that, in 
certain limits, the skew-scattering, contribution is sim- 
ply proportional to the resistivityEJIIiJ while the side- 
jump contribution is proportional to the square of the 
resistivity^. Then, we should have the simple expres- 
sion 

Ph = Pyx = ap xx + bp 2 xx , (2) 

which implies that the relative importance of these two 
contributions depends both on the temperature and on 
the impurity concentration. However, we show in this 
paper that, even if the relation (g) remains correct, 
the skew-scattering mechanism contributes also to the 
quadratic term in the case of impurity scattering. Such 
behaviour was already been shown by Kondorskii et alrn. 

The traditional way to calculate the anomalous Hall 
resistivity is to include the contribution of spin-orbit cou- 
pling in the transition probability (it leads to the skew- 
scattering provided one goes beyond the Born approx- 
imation) and in the velocity (it leads to the so-called 
anomalous velocity which gives the side-jump). While 
the skew-scattering can be obtained in a classical ap- 
proach it is claimed that the side-jump is a pure quan- 
tum effect. We shall discuss this point in the Sec. II 
of this paper. Most of the calculations of the anomalous 
Hall resistivity are based on the Bolztmann equation and 
used severe approximations, in particular cepcerning the 
side-jump contribution. Some calculations^ are based 
on the Kubo formalism, but surprisingly it is claimed 
that the side-jump contribution vanishes, and only the 
skew-scattering contribution is calculated. 

Although the anomalous Hall effect is an old phenom- 
ena which has motivated a lot of experimental and the- 
oretical studies, a unified model, able to calculate the 
skew-scattering and side-jump contributions on the same 
footing, was still missing. In this paper, we propose such 
a model. It is based on the Kubo formalism and has the 
peculiarity to be built from the Dirac equation. The jus- 
tification for such an approach is given in Sec. Ill where 
we discuss in detail two different approaches for solv- 
ing the anomalous Hall effect i.e., based on Dirac and 
Pauli equations, and study the consistency in the weak- 
relativistic limit of the expressions of the conductivity 
tensors obtained in these two approaches. In Sec. IV, we 
calculate the anomalous Hall conductivity of a disordered 
ferromagnetic bulk compound. The results are discussed 
in Sec. V. 



II. COMMENTS ON THE PHYSICAL NATURE 
OF THE SIDE-JUMP MECHANISM 

It is often believed that the side-jump is_a pure quan- 
tum effect and has no classical equivalents 2 ]. The usual 
description of the side-jump is then based on a quantum 
picture (see Fig. 0(b)) of a plane-wave transformed by 
scattering in the presence of spin-orbit coupling into a 
spherical wave whose center is shifted in a lateral direc- 
tion (perpendicular to the momentum and to the spin). 
The sign of the displacement is opposite for spin up 
(.s = 1) and spin down (s = — 1). A simple calculation 
in terms of phase-shift allows to determine this displace- 
ment. We start from the Pauli Hamiltonian 

2 

H = ^-- f i B {<T-B eff ) + W = H + W, (3) 

where cr is the Pauli matrix, TS e ff the effective magnetic 
field due to exchange interactions and W the total po- 
tential including the spin-orbit coupling 

W = V+-^(rTxVV)-p. (4) 

The state of the system \^>ks) after scatter- 
ing is given in the Born approximation by the 
Lippmann-Schwinger equation l^ks) = \k,s) + 
£ fc , s , W, s')G (k', s', e|)(V, s'\W\k, s), where e% and G 
are respectively the eigenvalues and the Green's function 
associated with Hq. The matrix elements of the potential 
are 

(*', s'\W\k, s) = VW (6 SS , + -^^(vs-s x k') • k) , (5) 

where Vkk' is the Fourier transform of V . As the spin- 
orbit term is imaginary, it will influence the phase of the 
spherical wave. Thus, for small spin-orbit coupling, the 
wave function ^E'fe s (r) = (r\^ks) which describes the wave 
after scattering can be expressed as 

* fcs (r) oc e lr k + ^^G (k' S ',4)y kk ,e^ k ', (6) 

k' s' 

where we have assumed that the effective magnetic field 
is along the z-direction. The center of the wave packet 
after scattering is given by 

h 2 

which is clearly spin dependent and means that the shift 
of the center of the wave packet is different for spin up 
and spin down. The lateral displacement, defined as 
5 s = r' s — r, is then equal to 

S s = ^(<r ss xp), (8) 

where we have introduced the length A which corresponds 
to the Compton wave length A c = fi/fuf in the case of 
free electrons. In real materials, Bergenia has shown that 
the spin-orbit coupling (i.e., A 2 ) is renormalized by band 
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structure effects by a factor a ~ 10 4 . We then obtain a 
lateral displacement S which is independent of disorder 
and of order A 2 fc^/4 « rs f0~ n m, in agree- 

ment with experimental results. Identical expression for 
6 was originally derived by Lyo et alB3. 

In a pure classical picture, such a lateral displacement 
can be experienced by a particle with spin. As a simple 
example, consider an electron with a charge e (e < 0) 
subject to a uniform electric field E = Eu x (E > 0) in 
the region x > 0; there is no field in the region x < 0. 
An incident electron coming from the region x < is 
reflected by the field as sketched in Fig. 0. 
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FIG. 2. Classical picture of the side-jump mechanism. The 
dashed line correspond to the non-relativistic trajectory of the 
particle and the solids lines to the relativistic trajectories for 
spin up (0) and spin down ((g)). 



The velocity is given by 
dH p 



v = 



dp 



Am 2 a 



■(<r x E), 



(9) 



and therefore contains an anomalous contribution v a = 
—eh(cr x E)/4m 2 c 2 arising from the spin-orbit interac- 
tion. In the field region (x > 0), where the trajectory is 
parabolic, the electron (we assume the spin to be along 
the z-axis) has an anomalous velocity along the y-axis, 
v v a = —eh((7 z E)/Am 2 c 2 . The electron therefore emerges 
with a shift along y, proportional to its spin a z . For an 
arbitrary electric field, the shift due to the anomalous 
velocity can be easily calculated 



S = 



v a dt = 



4to 2 c 2 



(erxE)cft, (10) 



with eEdt = dp, so that 
ha 



S = - 



Am 2 c 2 



dp 



^ X(P - P ' } ' 



(11) 



In the above derivation, we have assumed that the spin 
is perpendicular to the scattering plane. The lateral dis- 
placement that we obtain is consistent with the one ob- 
tained in the quantum picture. Indeed, the parallel can 
be simply done by replacing in this classical calculation 
the momentum by a momentum operator and by making 
the angular average over the final momentum p': thus 
(0) coincides with (|). 



III. COMPARISON OF THE DIRAC AND PAULI 
APPROACHES 

Generally, the calculations of the anomalous Hall con- 
ductivity are based on the Pauli Hamiltonian. However, 
in our modclization, i.e., within the framework of the 
Kubo formalism, it appears to be simpler to adopt a 
relativistic approach based on the Dirac equation. To 
justify that, let us first remember the derivation of the 
skew-scattering and the site-jump contributions in the 
Pauli approach. In presence of an exchange coupling, 
the Pauli Hamiltonian is H = H + H rc where H is the 
non-rclativistic Hamiltonian 



H 



2m 



[i B (<t ■ B eff ) + V, 



(12) 



and H rc the first relativistic corrections to the Hamilto- 
nian (order 1/c 2 ) 



H r 



P 



8m 3 c 2 Am 2 c 2 
% 2 

AV + H r 



(er x W) • p 



8m 2 c 2 



(13) 



which contains the relativistic mass correction, the spin- 
orbit coupling, the Darwin term and the relativistic cor- 
rection to the exchange coupling H rxc . Since the effect 
we are interested in results from the spin-orbit coupling, 
we do not need to give the explicit expression of H rxc 
(calculations and comments on this term are presented 
in Ref. |3l]). In this work, we do not consider the con- 
tribution of the periodic part of the spin-orbit coupling 
(i.e., due to the lattice) but only the aperiodic part due 
to the presence of impurities. In the Pauli approach, 
the velocity contains two parts. One resulting from the 
non-relativistic Hamiltonian v = p/m and another one 
resulting from the relativistic corrections 



P 2 P 

2m 3 c' 



h 



r(er x W) + v r 



(14) 



where v rac is the velocity related to H rxc . In this descrip- 
tion, the spin-orbit contribution to the velocity (second 
term in (Ji"4|)), the so-called anomalous velocity, appears 
in a natural and transparent way. When we insert this 
contribution in the Kubo formula, we obtain the side- 
jump contribution. It is also possible to isolate the spin- 
orbit contribution in the Green's function G associated 
with H by making the following expansion 



G — G + GH rc G + GH rc GH rc G 



(15) 



where G is the non-relativistic Green's function associ- 
ated with the non-relativistic Hamiltonian H. When we 
insert this expression in the Kubo formula and proceed 
beyond the Born approximation, we obtain the skew- 
scattering contribution. Therefore, in the Pauli approach 
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we get separately the skew-scattering and the side-jump 
contributions when the Green's functions and the veloci- 
ties are respectively corrected by the spin-orbit coupling. 

One important problem in the Pauli approach is to 
treat disorder. Actually, the spin-orbit coupling intro- 
duces two things: off-diagonal disorder (in the tight- 
binding approximation) and disorder in the velocity 
through the anomalous velocity. The second consequence 
is critical because it is then difficult to calculate pre- 
cisely the vertex corrections and accordingly the anoma- 
lous Hall resistivity. To avoid these problems, we have 
chosen to base our model upon the Dirac equation in- 
stead of the Pauli equatkww-,In presence of an exchange 
coupling, it has the formcSEa 



H = c(a-p) + j3mc 2 + V - [i B (i (tr ■ B eff ) 



(16) 



where the first term is the kinetic energy, the second term 
the mass energy, the third term is the potential and the 
last one the exchange coupling. From (16), we see that 
the velocity is simply 



dH 
dp 



er 
a 



(17) 



At this level, there appears an apparent contradiction be- 
tween the two approaches since, in the Dirac approach, 
contrary to the Pauli approach, we do not have any spin- 
orbit contribution to the velocity (anomalous velocity). 
It is therefore not clear a priori whether the side-jump 
mechanism would emerge from the Dirac approach. Ac- 
tually, in the Dirac approach, the spin-orbit coupling, 
although it does not appear explicitly, is properly taken 
into account. Therefore, the conductivity should con- 
tain simultaneously the skew-scattering and the side- 
jump contributions as well as higher order contributions 
in 1/c 2 . However, the expressions of the conductivity ob- 
tained in the Dirac and Pauli approaches should coincide 
in the weak-relativistic limit. To check this, we have cal- 
culated, in a formal manner, the weak-relativistic limit 



up to order 1/c 2 of the conductivity obtained from the 
Dirac equation and compared it with the conductivity 
obtained from the Pauli equation. The determination of 
the conductivity tensor is performed in the Kubo formal- 
ism. In certain limits, the conductivity can be expressed 
as a product of operators, namely Green's functions and 
velocities. However, the formulations propnsedJn the lit- 
erature are often confused or even wrongc3~t3 concern- 
ing the off-diagonal elements of the conductivity tensor 
due to an. abusive generalization of the Kubo-Greenwood 
formulae^. In order to clarify the situation, we present 
in Appendix A the derivation aLthe conductivity tensor 
from the original Kubo formulaH and summarize the dif- 
ferent stages apd approximations which lead first_tp the 
Bastin formulaO and finally to the Streda formulate. We 
show that the latter is a sum of two terms, dh and erf/ 
respectively given, in the limits of independent electrons 
appro ximat ion, zero temperature and zero frequency, by 
flAjj ) and ( |Ajj ) 



a 1 - 



„ - 47r -Tr( Vi (G+ - G-) Vj G- - ViG+v^G* 



Tr((G+-G-)(r i v j -r j v i )) c , 



-G-))e> 

(18) 



where i and j are the direction indices, f2 the volume of 
the sample, (...} c denotes the configurational average, G + 
and G~ are the retarded and advanced Green's functions 
at the Fermi level: G ± = G(e F ±i0) = (e F ± iO - H)^ 1 . 
The procedure that we follow is first to insert the Dirac 
velocity and Dirac Green's function in (^), next to per- 
form a weak-relativistic expansion of cry and finally to 
compare it with the expression obtained in the Pauli ap- 
proach. The Dirac velocity is given by (|i~7j ) and for the 
Dirac Green's function, we have used an exact expression 
derived from ([!(]) and given in Ref. |3l] by Eq. (A3) 



G ± = 



G i 



Zinc 



eff) 



' 2mc 



2mc 2 



(19) 



where the operators and are given by 



IP 



1 + 



LV £ F -V - H B (cr -B eff )) 



2mc 2 



1 



(rr-p)G ± ( ( T-p) 

2m 



(20) 



(21) 



The details of the calculations are presented in Appendix 
B. The determination of the conductivity is done up to 
order 1/c 2 . It is shown that the identification with the 
Pauli approach is successful only when one considers the 



total conductivity = aL + cr{J . Indeed, when we 
compare the expression of of - obtained in the Dirac ap- 



proach (sec (B5) for order l/c° and (B9) for order 1/c 2 ) 
to the one obtained in Pauli approach, we obtained dif- 
ferent terms which are exactly canceled by terms in a\j 



(see ® for order l/c° and (gig ) for order 1/c 2 ). The 
non-relativistic limit of the total conductivity obtained 
in the Dirac approach is 



_( ) = e a 



Tr 



^(G+ -G-p-G- 
m m 



4 



m m 



Tr((G+ 



G~ 



m m 



(22) 



which corresponds exactly to the conductivity obtained 
from (|is|) when one inserts the non-relativistic velocity 
v = p/m and the non-relativistic Green's function G. 
The last term in Eq. ( p2| ) is zero in absence of exter- 
nal magnetic field. The fact that a supplementary term 
in crjj ^ is present in the Dirac approch and not in the 
Pauli approach has serious consequences when one ne- 
glects ajj^ because it leads to an additional contribu- 
tion at order l/c° to the off-diagonal conductivity which 
does not disappear in the non-relativistic limit and thus 
would give unphysical results. At order 1/c 2 , the total 
conductivity obtained in the Dirac approach is 



~ (2) ~SS , ~SJ , -or 
a ij =(X ij + (T tj 



(23) 



where crf- s contains the terms which lead to the skew- 
scattering 

aff = -^Tr ( — {G + H rc G + - G~ H rc G~pG- 
J Anil \m m 

+— (G + - G~)—G~ H rc G~~ 
m m 

-PiG+H rc G + ^(G + -G-) 
in m 

G H TC G 



Pl G +Pl(G + H rc & 
m m 



(24) 



of/ contains the terms which lead to the side-jump 



e 2 h, 



G-p-G- 
m 

(v rc ),G+^(G+ -G-) + ^(& 
m m 



- ^G+(v rc ) J (G + 
m 

and af r is equal to 



G" 



-G-)(v re )jG- 
(25) 



-Tr ((G+-G-) {n{v rc ) 3 -r.j{v rc ) t ) 

G + H rc G+ - G-H rc G-\ (rP - r^)) . 

I \ m m/ I c 



(26) 



In addition to the skew-scattering and the side-jump con- 
tributions to the anomalous Hall effect, we identify a new 
contribution, 5°J , which is related to the orbital momen- 
tum L = rxp. The expression ( |23| ) of the conductivity 
corresponds exactly to the one which is obtained from 
( |l8| ) when one inserts the first order corrections to the 
velocity v rc and to the Green function GH rc G where v rc 
and H rc are given respectively by Eqs. ( |l3| ) and (14). 
We have then proved in the weak-relativistic limit (up 
to order 1/c 2 ) the coincidence of the conductivity in the 
two approaches. 



In summary, from the Pauli Hamiltonian, we get the 
skew-scattering and the side-jump contributions sepa- 
rately while, from the Dirac Hamiltonian, we get the both 
contributions and also higher order in 1/c 2 contributions 
simultaneously. Therefore, in a full relativistic Dirac de- 
scription, it will be difficult to assess the importance of 
each contributions. However, this approach has a great 
advantage over the Pauli approach: it allows a simpler 
treatment of the disorder because, in contrast to the Pauli 
approach where both the velocities and the Green's func- 
tions contain disorder, the disorder is only present in the 
Green's functions. It is thus possible to take one of the 
velocity operator outside of the configurational average 
and to calculate precisely the vertex corrections to the 
conductivity. For this reason, the Dirac approach should 
be more efficient to calculate the anomalous Hall resis- 
tivity. 

In the next section, we present a direct application of 
our model. In order to perform the analytical calcula- 
tions, we restrict ourselves to the weak-relativistic limit 
and to approximate calculations of the vertex corrections; 
then the results that we obtain can still be compared to 
the ones obtained from the Pauli approach. 



IV. ANOMALOUS HALL CONDUCTIVITY OF A 
FERROMAGNETIC COMPOUND 

In this section, we present the calculation of the 
anomalous Hall conductivity of a ferromagnetic bulk 
compound submitted to a potential. This calculation is 
done in both Dirac and Pauli approaches in order to show 
the similarities and the differences between these two ap- 
proaches. We consider a system with a cubic symmetry 
and a magnetization along the z-axis. Thus, the conduc- 
tivity tensor has the form 







@xy 
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a 



(27) 



We are only interested in the relativistic corrections to 
the off-diagonal elements which correspond to the anoma- 
lous Hall effect. We do not study the relativistic correc- 
tions to the diagonal elements which correspond to the 
anisotropic magneto-resistance (AMR) and lead to a dif- 
ference of order 1/c between a xx and a zz . Thus, in this 
work, the diagonal elements are calculated at order l/c° 
and by consequence are all equal, while the off-diagonal 
elements are calculated at order 1/c 2 . To get analytical 
expressions, we have made several approximations: free- 
electron approximation, weak-scattering limit and weak- 
relativistic limit for the Dirac approach. In Sec. Ill, we 
have shown that the conductivity is equal to 



- Vi G+v j (G + -G-)) c 



Tr( Vl (G+ -G-)v 3 G- 



5 



"4^ Tr < (G+ ~ G ^ inVj ~ ' (28) 

where the Green's function is associated with the to- 
tal Hamiltonian: G ± = (eFiiO-H)- 1 = (e F ±iO-H - 
W)~ l where Hq is the non-perturbed Hamiltonian and 
W the perturbation (equal to the potential V in the Dirac 
approach and to V + H so in the Pauli approach where 
H so is the spin-orbit coupling). The explicit form of the 
potential V does not enter in the calculations, thus the re- 
sults obtained below apply for both impurity scattering 
and phonon scattering in the adiabatic approximation. 
We modelize the compound in the following way: the to- 
tal volume of the sample fl — L 3 is divided into N cells of 
volume fio = a 3 . In each cell, the potential takes a con- 
stant value V with a probability distribution P(V) which 



is characterized by its moments (V n ) c = J P(V)V n dV. 
A proper choice of the energy origin yields {V) c — 0. We 
assume that there are no correlations in the value of the 
potential in different cells. In this first approach, we ne- 
glect in fl2S| ) the contribution of the terms which involves 
product of two advanced (or retarded) Green's functions. 
Such— an approximation is justified in the weak-disorder 
limitE3. In Appendix B, we have shown that afj, calcu- 
lated in the Dirac approach, contains two parts, the first 
one related to the orbital momentum, which is negligi- 
ble in our model, and the second one which is exactly 
compensated by terms in . Then, we do not need to 
calculate this contribution. The conductivity reduces to 

^ = ^^(^G + v 3 G-) c . (29) 




FIG. 3. Illustration of the conductivity with the help of Feynman diagrams. The total conductivity (hatched diagram), 
expressed like an infinite sum of diagrams involving the non- disordered Green's function Go (thin curve line), can be rewrite 
like an infinite sum of diagrams involving the average Green's function G (bold curve line). The wave lines refer to the velocity 
and the dashed lines to the potential. 



We introduce first the t-matrix T = W + WGqT which 
allows to write the Green's function as G = Gq + GqTGq 
where Go is the non-perturbed Green's function. Insert- 
ing this in (29), we get 



-Tr(v l G+TG+v J G TG-) c 



13 2nn 

e 2 h , 
+ 2tt£T 



(30) 



This equation can be illustrated with the help of Feyn- 
man diagrams as is done in Fig. [| The conductivity 
Oij , represented by the full diagram, is then express as 
a sum of an infinite number of diagrams. Only few of 
them are depicted in Fig. |^: diagram (a) which corre- 
sponds to the first term in ( |30| ) and diagrams from (b) 
to (f ) which are some representative samples of the kind 
of diagrams which give the second term in (^). The 
main approximation done in our calculation is to ne- 
glect the crossed diagrams which correspond to weak- 
localization corrections, i.e., we neglect diagrams such 
(e) and (f) and we keep only the so-called ladder dia- 



grams. Weak-localization corrections to the anomalous 
Hall conductivity are discussed in a separate paperEa. 
We introduce the configurational average Green's func- 
tion G = (G) c which can be written with the help of the 
self-energy E = (WG W) C + {WG a WG W) c + ... since 
G = (ef — Hq — S) _1 . When we neglect the crossed 
diagrams, (|3^) can be written as 



2ttH 

e 



2 h 



, )?rQ -Tr (v l G+T'G+v J G-T'G-) c , (31) 

with V solution of T = W+WGT'. The first term in the 
right side hand is the so-called bubble term (= d\^ bhle ) 
and the second one corresponds to the vertex corrections 
(= a^ rtex ). Within this transformation, the calculation 
of the conductivity is then reduced to two distinct prob- 
lems: determination of the average Green's function (i.e., 
the self-energy) and calculation of the vertex corrections. 
Because of the weak-scattering limit, we keep in the self- 
energy and the t-matrix the lowest sufficient orders 



G 



£ = (WG W) C 
V = W + WGW 



(32) 



In the t-matrix, we have to keep the terms up to the sec- 
ond order with V because it is necessary to go bewand the 
Born approximation to get the skew-scatteringU. The 
explicit calculation of j3l| ) in the approximations j3^ ) 
for the Dirac and Pauli approaches is presented in the 
next two sections. 



-bubble 



e 2 h 
2^0 



y~^(fc, s\vi\k, s')(k, s'\G + \k, s') 



kss' 



x (k, s'\vj\k, s)(k, s\G \k,s). 



(39) 



The configurational average (...) c has been dropped be- 
cause in the Dirac approach the velocity is a non- 
disordered quantity. At order l/c°, only the diagonal 
elements (s = s', no spin-flip) of the velocity and the 
particles in the upper band contribute, then we have 



A. Dirac approach 

We assume free electrons in a uniform effective mag- 
netic field B eJ ff parallel to the z-axis and submitted to a 
potential. The non-perturbed part of the Hamiltonian is 



H = c (a • p) + (J3 - l)mc 2 - fi B Po- z B 



effi 



(33) 



and the perturbation part is simply the potential W = V . 
The matrix elements of the average Green's function are 
{k,s\G^\k,s) = {ep — Ef. ± iTi/2r k )~ 1 where the eigen- 
values Ef, of ([33]) are in the weak-relativistic limit equals 
to 



2 1,2 



h'k 
2m 



1 



for the upper band, and 



4 = - 2mc +°l^)> 



(34) 



(35) 



for the lower band. The s index refers to the spin (s = 1 
for spin up and s = — 1 for spin down) , the k index refers 
to the upper band and the k to the lower band. The life- 
time t| which appears in the expression of the average 
Green's function is given by 



— = -Im(fc,s|S + |/c,s) = -lm(k,s\(VG£V} c \k,s) 



2r, 



= irn Af s (e 8 k ){V 2 ) c , 



(36) 



where J\ s is the density of states of spin s by unit volume. 
In the Dirac approach, the velocity v is simply equal to 
col. Because we have chosen to work in the basis where 
the non-perturbed Hamiltonian Hq (and by consequence 
the Green's function Go) is diagonal, we have to calculate 
the velocity in this basis, we get 



u(k) 
ccr - 



ccr + o[ 

o(c°) 



where u(k) is the (2 x 2) matrix 



Kk 

u(k = — 

m 



1 
1 



(37) 



(38) 



We have now all the ingredients to calculate the bubble 
term in (till) 



-bubble 



e 2 h 3 



(40) 



The dispersion law e 8 k given by (gj) is isotropic at order 
1/c . Then, the angular dependence is entirely contained 
in the factor kikj, which means that only diagonal com- 
ponents of the conductivity are different from zero. To 
order l/c°, the vertex corrections to the diagonal compo- 
nents vanish, so that the total conductivity da is equal 



to a\f hble . After integration over k, we get 



2 Afi 



(41) 



which corresponds to the Einstein relation with two spin 
channels where Z s = v f t f is the mean-free-path, v F and 
Af s are the velocity and the density of states by unit vol- 
ume at the Fermi energy for spin s respectively (identical 
expressions are obtained for a yy and a zz ). The diagram 
which gives this contribution is depicted on figure |4| 




FIG. 4. Bubble diagram contributing to the diagonal con- 
ductivity a xx . The signs +/— refer to the retarded/advanced 
average Green's function G . 

The off-diagonal components of the conductivity arise 
only when we take the vertex corrections into account. 
If we expend the t-matrix up to the second order in V, 
we get 



from 



-vertex 
ij 



-TT(viG + {V + VG + V)Ct 



2irfl 

xvjGr(v + vcrv)G-) c . 

We then need the potential in the new basis 



(42) 



o(c°) 



(43) 



where U(k, k') is the (2 x 2) matrix 
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E/(k,k') 



1 - 



ri 2 ((k'-k) 2 +2i(k'xk)-e z ) 



^((k^+ik ')k' z 



8m 2 c 2 

-(kj;-\-iky )k z — z(k' xk) ■ (e K +2e y )) 



4m 2 c 2 



fi 2 {(k' so — ik' y )k' r — (k !c —ik y )k z — i(k' xk)-(e.,. —ie y )) 
4m 2 c 2 

-, _ ft 2 ((k'-k) 2 -2i(k'xk)-e^) 
1 



(44) 



and V(q) — J dr e tcl ' r V(r)/Q is the Fourier transform of 
the potential. When we study in detail all the diagrams, 
included in (|4^) , we see that only two kind of diagramsE-3 
contribute to the conductivity at order 1/c 2 . These dia- 
grams are depicted on figures and ^ (left column). 





FIG. 5. Diagrams contributing to the off-diagonal conduc- 
tivity a xy through the skew-scattering mechanism in both 
Pauli and Dirac approaches. The number in bracket indicates 
the order with 1/c of the matrix elements of the velocity (wave 
line), average Green's function (bold curve line) and potential 
(dashed line). It is omitted in case of zero order with 1/c. 

The first series of diagrams (see Fig. |^) involves ve- 
locities at order l/c°, Green's functions at order l/c°, 
which means that only particles in the upper band con- 
tribute, and potential twice at order 1/c and once at 
order 1/c 2 which ensures a total order of 1/c 2 for the 
conductivity. The diagrams of this first series correspond 
to the skew-scattering mechanism. The second series of 
diagrams (left column in Fig. |^) involves one velocity at 
order l/c° and one at order c, three Green's functions at 
order l/c° and one at order 1/c 2 , which means that we 
have a transition of particles between the upper and lower 
bands, and potential one time at order l/c° and one time 
at order 1/c which ensures a total order of 1/c 2 for the 
conductivity. The diagrams of this second series corre- 
spond to the side-jump mechanism. In the following, we 
present the explicit calculation of the conductivity due to 
these two series. Let us start with the skew-scattering. 
We present the calculation of the diagram (a) in Fig. 
which gives 



a 



e 2 h 



{(k,s\v x \k,s)^(k,s\G + \k,s) 



(0) 



x (k, s\V\k', s) (0) (k r , s\G + \k', s) (0) (k', s\V\k", s) (0) 

x (k", s\G + \k", s) (0) (k", s\v y \k", s) (0) (&", s\CT\k", s) (0) 

x (k",s\V\k, S )W(k,s\G-\k,s}(°A . (45) 



The number in bracket indicates the order with respect 
to 1/c of the matrix elements like in Fig. || when the or- 
der is different from zero. We remark that for a total 
order 1/c 2 of the conductivity, the spin is conserved dur- 
ing the process (no spin-flip scattering). We insert in this 
expression, the matrix elements given by (pq) and (E3) 



and perform the integration over fc, k' and k" . The final 
contribution to the conductivity corresponding to the di- 
agram (a) is a complex quantity. The calculation of the 
diagram (b) gives the conjugated expression, then the to- 
tal contribution due to the skew-scattering mechanism is 
a real quantity equal to 



~ss 



ttto 2 A 2 (V 3 ) c 

6h 2 (v 2 ) c 



(•A/- T flofrL(4) 



~SS1 
"xy 



xy 



(46) 





VA/ (c) 

'(2) y 



~aJK 2) (-b ... 

w (d) 



FIG. 6. Diagrams contributing to the off-diagonal conduc- 
tivity o xy through the side-jump mechanism in Dirac ap- 
proach (left column) and Pauli approach (right column). The 
number in bracket indicates the order with 1/c of the matrix 
elements of the velocity (wave line), average Green's function 
(bold curve line) and potential (dashed line). 

We turn now our attention to the side-jump mecha- 
nism. The diagram (a) of Fig. ^| gives 



u xy 



e 2 h 



Q (< fc > s\v x \k, -s)^{k, s\G+\k, -s) 

kk' s 

x (fc, -s\V\k', s)« (k', s\G+\k', «}(°) {k\ s\v y \k', s) (0) 
x(k', s\GT\k', s) (0) (fc', s\V\k, s) (0) (fc, s\G~\k, s) (0) 



(2) 



(47) 

In this mechanism, due to the presence of off-diagonal 



elements in the velocity (|37|) and the potential 
particle of the upper band e s k experiences a virtual tran- 
sition in the lower band e^ s associated to the opposite 
spin. We perform the integrations over k and kl , add 
the contributions of the four diagrams (^|a) to (^|d) and 
finally obtained 
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~SJ 2>r 2S ^i 2.r 2Siv F 

a s J = —e M\ — + e Mi — 



'xy 



= a SJT 



u xy ' 



(48) 



where 6 s is the transverse displacement (or side-jump) 
given by hvp/Amc 2 = X 2 k F /4. The expression of a xy for 

) of d xx but 
In contrast 



the side-jump is similar to the expression m 
instead of the mean-free-path I, we have 2t 



to a X y , the side-jump contribution to the off-diagonal 



conductivity is independent of disorder. 

In the case of a parabolic band, the Einstein relation 
( ffl| ) reduces to the Drude formula with two spin channels 



to 



i n l' 



(49) 



where n s — raJ\f s (v F ) 2 /3 is the electron density for spin 
s. The skew-scattering 
tions yield 



3) and side-jump ([48|) contribu- 



~ss 

xy 



e 2 A 2 ^ (^ 3 ) c 

" 2h h(v 2 ) c 



(n 2 f, 



and, 



~SJ 
xy 



2 \2 



e 2 \ 
~2lf 



■(n T -raj.). 



(50) 



(51) 



B. Pauli approach 

In the Pauli approach, the Hamiltonian given by (||) is 
the sum of a non-perturbed part and a perturbation W 
given by (|j) which contains the potential and the spin- 
orbit coupling. The velocity associated with this Hamil- 
tonicn consists of a normal part and an anomalous part 
due to the spin-orbit coupling 



and gives the same final expression ( }46| ) as the weak- 
relativistic limit of the Dirac approach. 

Concerning the side-jump, the correspondence between 
the two approaches is not so simple. In the Dirac ap- 
proach, we have seen that a virtual transition occurs from 
the upper band to the lower band. In the Pauli approach, 
no such a transition can take place because there is only 
one band. However, we have a supplementary part in 
the velocity, the anomalous velocity which is of order 
1/c 2 and leads to the side-jump mechanism. The corre- 
sponding diagrams are depicted on the right column of 
Fig. ^|. For each diagram in the left column (i.e., in the 
Dirac approach), we have an equivalent diagram in the 
right column (i.e., in the Pauli approach). The change 
between the left and right column corresponds to a vertex 
renormalization because the matrix elements of the prod- 
uct vGV in the Dirac approach are equal to the matrix 
elements of v in the Pauli approach 



Kk 

(k, s\v\k',s') = — 8 kk '5 ss i 
m 

+ V(k' -k)- 



ih 



s , x (k - k') . (53) 
Thus, when we calculate, for example the diagram (^a') 
ff 8° = TTi £ ((k, S \vx\k',s)^(k', S \G+\k' >S )^ 



kk' s 



x(fc',s|i; y |fc',s} (0) (fc',s|G-|fc',s) 



(0) 



x (k',s\V\k,s) {0) (k,s\Gr\k,s) 



\(0) 



(54) 



we obtain the same contribution than from the expres- 
sion ( p!7| ) of the diagram (||a). The final result, after 
summation over the four diagrams (^|a') to (^|d'), is then 
identical to (El) . 



to 4m 2 c 2 



(<T X W) 



(52) 



We neglect the contribution of the spin-orbit coupling 
in the life-time, then the average Green's function is 
G ± = (sp — H o =F ifi/r^ 1 where is given by (36). 



Because of this approximation, the derivation of the di- 
agonal conductivity a xx is similar to the one done in the 
Dirac approach and we obtain the expression ( f4l| ) . 

The off-diagonal elements of the conductivity are ob- 
tained from the vertex corrections. For the skew- 
scattering, the diagrams which contribute are exactly the 
same than in the Dirac approach (see Figj5) because the 
only matrix elements of the potential (J44) which con- 
tribute in the Dirac approach correspond precisely to the 
matrix elements of the potential in the Pauli approach 
(see Eq. (|j)). The other terms in ( [44] ) are Darwin- like 
terms and do not contribute to the off-diagonal conduc- 
tivity. Then, in the Pauli approach, the skew-scattering 
mechanism corresponds to the same Feynman diagrams 



V. DISCUSSION 

We now briefly discuss the influence of impurity scat- 
tering and phonon scattering on the resistivity and on 
the anomalous Hall resistivity, which are, in the limit 
o xy < cfxx simply given by p xx ~ l/a xx and p H = 
—p xy ~ &xy/d' xx - The only terms which depend on the 
scattering in the expressions of a xy given by (Eq) and 



(g§) and of a xx given by (g 
and (V" 3 ) c . Indeed, we have 



]), are the moments (V ) c 



xy 



X 



(55) 



indep. of {V n ) c 



Then, the variations with the moments (V 2 ) c and (V 3 ) c 
of the resistivity and the anomalous Hall resistivity are 
like p xx cx (V 2 ) c , pH ex (V 3 ) c and pgj ex (V 2 ) 2 C . 
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To illustrate the dependence with disorder in the case 
of impurity scattering, we consider a binary alloy A X B\- X 
for which 



and, 



{V 2 ) c = x(l-x){e A -e B f 



{V 3 ) c = x(l-x){l-2x)(e A ~e B ) 3 , 



(56) 



(57) 



where £a(B) is the value of the potential V on site A{B) 
and x is the concentration of sites A. Keeping the lowest 
orders in x (weak-disorder limit), we get 



Pxx X 

Pxy <X(X- 3X 2 ) 
Pxy X 2 



(58) 



which is in agreement with the simple relation given by 
(||) but in contradiction with the common belief that 
the quadratic term would arise only from the side-jump 
mechanism. In fact, the skew-scattering mechanism, 
which is responsible for the linear term, gives also an 
important contributiop_to the quadratic term, a result 
that Kondorskii et alr3 have already obtained. In ad- 
dition, our calculations specify all the approximations 
which founded the relation (^) and show that it should 
not be valid in the general case, in particular for high- 
disordered system, high-relativistic limit, complex band 
structure or in the case of heterogeneous systems, such 
as thin films or multilayers. 

Something may also be said about phonon scattering. 
Due to the fluctuating sign of the potential generated by 
phonons, the third moment {V)l can be expected to be 
very small and accordingly, the skew-scattering contri- 
bution (^6|) to the conductivity is negligibleE 2 ]. We have 
then 



Pxx oc (V 2 ) c 
P^J - o 
PU oc (V 2 )l 



(59) 



which yields the simple relation p xy oc p xx , in agreement 
with experimental resultsEJ. 

The Hall angle, which corresponds to the angle be- 
tween the electric field and the charge current, is an 
important quantity. For an applied electric field in 
the x-direction and an effective magnetic field in the z- 
direction, we have tg(0n) = j y /j x = o yx /a xx . The con- 
ductivity elements a xx and a yx are in a first approxima- 
tion the sums of contributions due to spins up and down. 
We can thus define a spin-dependent Hall angle 



tg(0 



-•TU) 
Jy 

,•1(1) 

Jx 



u XX 



(60) 



We insert the expressions ([l6]) and ( |48| ) of the skew- 
scattering and side-jump off-diagonal conductivities as 
well as the expression ([II]) of the diagonal conductivity 
and obtain for spin s 



'H 



25 s 



irm 2 \ 2 (V 3 ) c 
6h 2 (V 2 ) t 



(61) 



^ are not only opposite in sign, they take dis- 



6 H and w H 

tinct absolute values due to spin polarization. As a conse- 
quence, the spin current (j^ — j^) has a longitudinal (i.e., 
along the x-axis) and a transverse (i.e., along the y-axis) 



component and the charge current (j^ + j^) acquires a 
transverse component which corresponds to the anoma- 
lous Hall effect. In a paramagnetic material, Eq. ( pl| ) 

T I — 

yields H — —6 H and both the transverse component of 
the charge current and the longitudinal component of the 
spin current vanish. However, the transverse component 
of the spin current remains. It corresponds precisely to 
the "spin Hall effect" recently proposecO. 

In the case of impurity scattering and in the weak- 
disorder limit, the magnitude of the Hall angle is deter- 
mined mostly by the skew-scattering contribution. In- 
deed, in this limit we have H J w 28/1 w 10~ 3 rad 
whereas Off « 7r(l - 2x)(ea — sb)sf /SnuPW ~ 
5.10~ 2 rad where we have taken I w 200 A, x m 0.2, 
ea ~ £b ~ 2 eV, £f ~ 10 eV, the band width W w 5 eV, 
mc 2 ~ 500 keV and a band factor a « 10 4 . For simplic- 
ity, we have dropped the spin index. This ordex.of mag- 
nitude is consistent with experimental resultsEHl. When 
the disorder increases, the mean-free-path I decreases sig- 
nificantly which means, since the quantity S is disorder 
independent, an increase in the side-jump contribution 
to the Hall angle. However, the skew-scattering contri- 
bution to the Hall angle increases in the same way. It is 
thus not possible to predict in this first approach which 
contribution is dominant in the high-disorder regime. 

In the case of phonon scattering, the Hall angle con- 
tains mostly the side-jump contribution Oh = H J fa 25/1 
which is of order w 10~ 2 rad where we have used S « 
10~ n m and I = tvf with a relaxation time r « 10 -15 s 
and vf ~ 10 6 m.s -1 . 

To summarize, we have, in this article, proposed a 
model based on the Dirac equation and on the Kubo for- 
malism which allows to calculate on the same footing the 
anomalous Hall conductivity due to both skew-scattering 
and side-jump mechanisms. The consistency of this ap- 
proach with the one based on the Pauli equation has 
been studied in detail in the weak-relativistic limit. In 
particular, we have shown that in order to calculate the 
anomalous Hall conductivity one has to consider in the 
Dirac approach the total conductivity 5^ + a(j , other- 
wise unphysical results are obtained. Next, we applied 
our model to treat a disordered ferromagnetic bulk com- 
pound submitted to a potential in the free electron ap- 
proximation, weak-scattering and weak-relativistic lim- 
its. By these means, we have obtained explicit expres- 
sions for the anomalous Hall conductivity for both skew- 
scattering and side-jump mechanisms (given by (fl6l ) and 
(ft8|))- In addition, we have highlighted the difference con- 
cerning the Feynman diagrams describing the side-jump 
mechanism in the Dirac and Pauli approaches and have 
shown that it corresponds to different vertex renormal- 
izations. 
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APPENDIX A: FROM THE KUBO FORMULA 
TO THE STREDA FORMULA 

In the linear response approximation, Kubo has shown 
that the conductivity, tensor is related to a two currents 
correlation functionE3 



dij{ui) = f2 lim 







+ 00 



xTr( poJj^Jiit + ihX) 



where it is assumed that the applied field leads to 
a time-dependent perturbation of the form: H'(t) = 
-ffpexp(^(— ?kj + is)), n is the volume of the sample, 



(3 = po is the density matrix in equilibrium in 

absence of perturbation, Ji is the i-component of the cur- 
rent density operator in the Heisenberg representation 
and (...} c denotes the configurational average. Follow- 
ing LuttingercS, we obtain in the independent electrons 
approximation 



_ A(t+ih\)(s n -e m ) 



(n\Ji(t + ih\)\m) = e 
(Al) where we have used H = ^ n e n c 



n\Ji\m), (A2) 



and defined 



J as the current density operator in the Schrodinger 
representation. Using the relation Tr[poOm a n a p a q ~\ = 
5 mq 5 np f(e m )(l - /(£„)) where /(e) is the Fermi-Dirac 
distribution function, we get 



? ii (w) = n _lim / d\e- x ^~^ dtJ2{f(e m )(l- f(e n ))e^ 



.0+ 



^ ( — 7w*>+ is-\-E n — e r , 



m\Jj\n){n\Ji\m)j . (A3) 

The integration over A leads to a factor (1 — e _ ^ e "~ £m ))/(£„ — e m ) which can be simplified with f(e m ) (1 — f(s n )) as 

(A4) 



-/(£m) (1 - f(Sn)) = 



Inserting this in ( |A3| ) and performing the integration over t, we obtain 

/(em) - f(£n) 



aa(uj) — ihVl lim > 



(£n - £m)(£n - £m ~ flW + is) 
I 



m\Jj\n)(n\Ji\m) 



(A5) 



We shall now make some transformations of this expres- 
sion in order to get the Bastin formula. 

We restrict our derivation to zero frequency (from now, 
we drop the ui variable). After inserting the identity 
J+°° deS(e-H) = l in @, we obtain 



/+00 
de 

-°° nm 

f(e)S(e-e n ) 
(e - e m )(e - e m + is) 

We remark that 

1 



f(e)S(e 



(e„ - e)(e n - e + is) 
m\Jj\n)(n\Ji\m) \ . (A6) 



lim 



lim — 



^0+ (e n — e)(e n — e + is) s^o+ de \e n — e + is 



(A7) 



then we have 
a a = —ihQ lim 



s-»0+ 



+ 00 



de \ e 



def(e) 
1 



(n\Ji\m)S(e 



-{m\Jj\n)S(e - e n ){n\ Ji\m) — 



de V £ 



which can be expressed as 



n 



<1:J cm {v^^-VjSiB-H) 



ViS(e — H)v 



dG~{e) 
de 



(A8) 



(A9) 



where we have introduced the Green's function G ± (e) — 
lim s ^ + (e — H ± is)^ 1 and the velocity through the re- 
lation J = —ev/Q. This expressioiufor the conductivity 
was first obtained by Bastin et alEU but in the particu- 
lar case of a Schrodinger Hamiltonian and made explicit 
use of the form taken by the velocity operator in the 
Schrodinger case. The present derivation is more general 
in the sense that it is independent of the explicit form 
of the velocity operator and is therefore valid both for 
the Schrodinger, Pauli and Dirac cases. The only restric- 
tion is the independent electrons approximation. This 
formula, called Bastin formula, is interesting because it 
expresses the conductivity as a product of velocities and 
Green's functions. However, it is still difficult to calculate 
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because of the integration over the energy e. By making 
an integration by parts, a factor df(e)/de appears instead 
of the factor f{e) and the integration interval will be thus 
reduced. 



In (A£), we express the delta function in terms 



of Green's functions using 5(e — H) = — (G + (e) — 
G~(s))/2iir. We keep one half of this expression and 
make an integration by parts on the second half, then we 
get 



e 2 h 



+ 00 



de 



dm 

de 



Tr(vi (G+(e) - G~{e)) VjG-(e) - ViG+^vj (G+(e) - G~(ej) ^ 



+ 



e 2 H 
4tt^ 



de 



de 



de 



(A10) 



The second term in this expression can be simplified 
by using the relations dG ± (e)/de = — (G ± (e)) 2 and 
ihVi — [ri, H] = —[ri, G _1 ] and by performing once more 
an integration by parts. Finally, the conductivity can be 

07/ where 



written as a sum of two terms = afj 



4 



e 2 h 
~4v^ 



—00 



de 



xTr/vi (G+(e) - G~(e)) 
-v l G+(e)v 1 (G+(e)~G-(e))\ , 

I c 



(All) 



and 



ij 



xTr 



+ - d£ dm 

-no de 



((G+(e)-G-(e)) (r^-r^)) • 



(A12) 



( |All| ) .and ( [A 12] ) correspond to the formula obtained by 
StredaES in the Schrodinger case. The present derivation 
shows that it holds also in the Pauli and Dirac cases. 
For the diagonal components of the conductivity tensor, 
&U is eoiial to zero and we obtain the Kubo-Greenwood 
formulae*] 



de 



dm 

de 



L (G+(e)-G-(e)) 
xv i (G+(e)-G-(e))^ . 



(A13) 



At zero temperature, the factor df(e)/de is equal to 
—S(e — £p), only electrons at the Fermi level contribute 
to the conductivity (for both diagonal and off-diagonal 



components). In conclusion, at w — and T = 0, the 
conductivity tensor can be expressed as a sum of two 
terms an = afj + aU with 



° 13 ~ 4t70 



and, 



Tr(vi (G + -G-)vjG- 
ViG+vj (G+ - G~) \ 

I c 



Tr(^(G+-G-)(r i v j -r j v i )j 



(A14) 



(A15) 



where we have dropped the energy reference £f by in- 
troducing the Green's functions at the Fermi level G = 
G(e F ± iO) = (e F ±i0- H)~ l . 



APPENDIX B: STREDA FORMULA IN THE 
WEAK-RELATIVISTIC LIMIT 

In this appendix, we give the detail of the calculation 
concerning the weak-relativistic expansion of the Streda 
condu ctivity starting from the Dirac equation. From 



(A14), we see that a\j is a combination of terms such 



as 



e 2 h t 
4t7^ 



Tr( ViG(z\)v : jG(z 2 ) 



(Bl) 



where Z\ and z 2 are equals to £i?±iO. When we insert the 
Dirac velocity ( |l7| ) and the Dirac Green's function ( |l9| ) 
in (Bl), make the explicit product of the 4 operators and 
take the trace over the lower and upper components of 
the wave function, we obtain the general form: 



Ay(*i,*0 = ^Tr(a i D(z 1 )^G(z 1 )* j D(z 2 )^G(z 2 ) 
+a l G(z 1 )^Q-\z 1 )D(z 1 )Q(z 1 )a 3 G(z 2 )^Q- 1 (z 2 )D(z 2 )Q(z 2 ) 
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+a i D(z 1 ) 



(B2) 



Similarly, when we insert ( p7| ) and ( |l9| ) in ( A15| ), we get 



4inn \\ 2m 



(Q+)- i ^+Q+ + D+^G+ - G-^-{Q-)- l D-Q- - D 



2m 



2m 



<* P ^ 

2m 



(B3) 



Expressions (B2) and flB3| ) are exact expressions with- 
out any assumption on the value of c. We will now calcu- 
late the weak-relativistic expansion of these expressions 
at orders 1/c and 1/c 2 in order to compare them with 
the expression obtained from the Pauli approach. 



in the Pauli approach and should not appear when we 
take the non-relativistic limit in the Dirac approach. In 
fact, we show below that this term is exactly cancelled 

Replacing D by 1 in (jB§), 



by an opposite term in cf*- 



we get of 7 at order l/c° 



1. Dirac conductivity at order l/c° 

In the non-relativistic limi t, D is simply equal to the 
unit matrix (see (|2C|)), thus (B2) can be rewritten as 



^Tr(^G( 
AttVL \m 



:ij^G(z 2 ) 
m 



+ -j— \ ai<TjG(z 2 ) + aj<7iG(zi) 



(B4) 



a 11 ^ 



-Tr 



G^-G-)(n p -i-r^ 

/ V m m 

e 2 h 



Eijk 



-Tr 



(<*( 



G + -G" 



(B6) 



where we have used the relations (or ■ A) (er • B) = 
(A ■ B) + icr ■ (A x B) and [r^/jj] = ihSij. The second 
term in the right hand side cancels the supplementary 
term in (B5) and we obtain finally for the total Dirac 
conductivity at order l/c° 



where we have used the fact that <7j (<x • p ) + ( tr ■ p) at — 
2pi. When we insert this expression in (A14), the con- 
ductivity o\a at order l/c° is 



_/(o) _ e n 
° l 3 ~ 4tt^ 



Tr 



Ei(G+-G-)^G- 
m m 



_P±Q+P]_(ri+ 



+£ 



i-jk 



4AirmQ 



Tr ( a k (G+ - G~ 



(B5) 



where Sijk = 1 if k} = {x,y,z\ or cyclic permuta- 
tions and Eijk = otherwise. This factor is introduce 
through the term (pi<jj — ctjOi) = 2i£ijkOk- The first 
term on the right hand side corresponds exactly to the 
contribution that we get in a non-relativistic description 
because in this case the velocity is v = p/m and the 
Green's function is simply the non-relativistic Green's 
function G. In contrast, the second term is not present 



-(0) ~/(o) 
a • ■ = (T ■ 

w %3 



-^-Tr( — (G^ 
47rO \ m v 



4i7rSl 



-Tr 



(g+ 



G')^G~ 



Pj 



m 

Pi 



— (G -1 
m 



G~ 



G~ ) in— —n— 



(B7) 



which corresponds exactly to the total conductivity 
obtained from (A14) and (A15) when we insert the 
non-relativistic Pauli velocity v = p/m and the non- 
relativistic Pauli Green's function G. 



2. Dirac conductivity at order 1/c 2 

To get &ij at order 1/c 2 , it is necessary to take into 
account the next terms in the expansion of D given by 
&■ D(z) » l-Q(z)(z-V-y B ((T-B ef f))/2mc 2 . Thus, 
from (JB2J), we get 



e 2 h , 
4^V 



-Tr 



-G( Zl )(v rc )jG{ Z2 ) + (v rc ) J G(z 1 )^G(z 2 ) 
m m 



+ Pi G (z 1 )^G(z 2 ) J ff rc G(z 2 ) + ^Giz^HrcGiz^Gizz) 
mm m m 



8m 3 



^cr i cr :) G(z 2 )cr ■ p(z 2 - V - /i B (er • B eJ y))<T • pG(z 2 ) + (Tj(TiG( Z \)a- ■ p(zi - V - ^sO ■ B eff ))cr ■ pG(zi) 
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~ ' 2m 3 c 2 M3 + < ^( 2:2 )) _ 4^2 ( CT '( Zl ~ F _ Vb{(t ■ B eff ))a j G{z 2 ) + (Tj(z 2 - V - fi B (cr ■ B eff ))cr i G(zi)) 

- 4m ' 3c 2 ( z i _ z 2) ((p x a)iG(zi)pjG(z2) -piG(zi)(p x ct^G^)) y , (B8) 

where H rc and v rc are the relativistic corrections at order 1 /c 2 to the Hamiltonian and the velocity respectively given 
by ( |l3| ) and (H)- The last term on the right side hand does not contribute because {z\ — z 2 ) — ► 0. Inserting the 
expression (|B8[) in the conductivity (A14), we get a-j at order 1/c 2 

4 (2) = ^Tr/(v rc ) i (G+ - G-)^G- - (v rc ) i G+^(G+ - G") + ^(G+ - G^v^C" 

--G + (v rc ) 7 (G+ - G~) + -{G+H rc G + - G-H rc G-)^G- + -(G+ - G~)—G~ H rc G~ 
m m m r- 

Pi r<+ tt rr+Pi rfi+ r<-\ Pi fi+Pi (n+ti n+ n-n r<-\ 

(_T il rc (_x l(_r — (_r I Lx (Lx il rc (_r — (_r tl rc \jc ) 

m to mm 

-e-tjfe i6j 7rTO 3 c 2^ Tr \ ( G+cr ' P( e -F - V - ^b(<t • B ej? ))cr • pG + - G a ■ p(e F -V- /jb(<x • B Gjff ))cr • pG J <j k 

+ i6Wc 2 n Tr ( (6+ " ^^ £F - y - ^ ' Be ^))^ ^ ff ^ £F ^ y ^ ^ • B «sff))°-o) • (B9) 

The first term corresponds exactly to the relativistic corrections that we get at oder 1/c 2 in the Pauli approch. The 
two last terms are supplementary terms which should not appear. We show that they are cancelled by terms in aj^ 2 \ 
Indeed, when we expand D up to the second order in 1/c in the expression (B3) of a(J , we obtain 

4 7(2) = -^ Tr ( (G + H rc G + - G~H rc G) (rg - ) + (G+ - G~ ) (n{v r c)j - rj(v r c)i) 



771 777 



+£ijk 1677rm 3 c 2 ^ Tl \ { G+(T ' p ( £F ~ v ~ ^b{it- B eff ))cr ■ pG + — G a ■ p{e F - V - /_t B (er • B eff ))(T pG ) <r k 

Tr/(G+ - G-){ai(e F -V- ti B (cr ■ B eff ))a 3 - Oj{e F - V - fi B (cr ■ B eff ))a z )\ . (BIO) 



16nm 2 c 2 Q \ 



The two last terms on the right hand side cancel the supplementary terms in ( B9 ) and we obtain finally for the total 
conductivity at order 1/c 2 

*$> = *T + ^ (2) = ^r ((v, c MG+ - G-f^G- - (v rc)i G+|(G+ - G-) 
+ ^1(G+ - G-)(v rc ),-G- - -G+(v rc ),(G+ - G-) + -{G+H rc G+ - G~ H rc G~)—G~ 

m m 777 777 

+ Pl(G+ - G-)^G-H rc G- - p ±G+H rc G + ^(G+ - G") - ^G + ^(G + iT rc G + - G~ H rc &- 
to m mm mm 

D 2 



' -Tr 



477rfi 



((<2 + H rc G+ - G- H rc G-) (rA-rj^j + (&+ - G") (r 4 (v rc ), - r^v^i)) • (BH) 



This expression corresponds exactly to the one that we 
get from the Pauli approach when we report the Pauli 
velocity at order 1/c 2 , v rc , and t he P auli Hamiltonian at 



order 1/c 2 , H rc , in (A14) and ( A15| ). For higher order 
terms (say of order l/c 2ra with n > 1), we can predict 
that similar cancellations occur when we consider the 
total conductivity a 



~I(2n) 



i J 



~II{2n) 

a v ■ 



We have thus proved 
in this appendix the consistence between the Pauli con- 
ductivity and the weak-relativistic limit of the Dirac 



conductivity. 

* Electronic address: crepieux@mpi-halle.mpg.de 
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